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ABSTRACT 



Context. Quasi-periodic propagating disturbances in coronal structures have been interpreted as slow magneto-acoustic 
waves and/or periodic upflows. It is important to understand the nature of these disturbances before proceeding to 
apply them. 

Aims. Here we aim to understand the nature of these disturbances from the observed properties using a three-hour 
imaging sequence from AIA/SDO in two different temperature channels. We also compare the characteristics with a 
simple wave model. 

Methods. We searched for propagating disturbances in open-loop structures at three different locations; a fan loop- 
structure off- limb, an on-disk plume- like structure and the plume/interplume regions in the north pole of the sun. In 
each of the subfield regions chosen to cover these structures, the time series at each pixel location was subjected to 
wavelet analysis to find the different periodicities. We then constructed powermaps in three different period ranges, 
short (2 min - 5 min), intermediate (5 min - 12 min), and long (12 min - 25 min), by adding the power in the individual 
range above the 99% significance level. We also constructed space-time maps for the on-disk plume structure to estimate 
the propagation speeds in different channels. 

Results. We find propagating disturbances in all three structures. Powermaps indicate that the power in the long-period 
range is significant up to comparatively longer distances along the loop than that in the shorter periods. This nature is 
observed in all three structures. A detailed analysis on the on-disk plume structure gives consistently higher propagation 
speeds in the 193 A channel and also reveals spatial damping along the loop. The amplitude and the damping length 
values are lower in hotter channels, indicating their acoustic dependence. 

Conclusions. These properties can be explained very well with a propagating slow- wave model. We suggest that these 
disturbances are more likely to be caused by propagating slow magneto-acoustic waves than by high-speed quasi-periodic 
upflows. We find that intensity oscillations in longer periods are omnipresent at larger heights even in active regions. 
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1. Introduction 

Small-amplitude propagating intensity disturbances 
are often observed with imaging telescopes in 
open/extend ed loop structures close to active regions 
on the sun dBerghmans fc Clettd [l99l iDe Moortel et all 
120001: iKirTg et all l2QQ3h and also in polar p lumes 
(iDeforest fc Gurmanl 119981 : iMdntosh etHI 120101 : 
Krishna Prasad et all 1201 lh . These are identified as 
alternating slanted ridges in the space- time maps, which 
are thought to be signatures of slow magneto-acoustic 
waves owing to some of the observed properties such as 
periodicities and propagation speeds (jDe Moortell 120091 
and references therein). Assuming them to be slow- mode 
waves, their contribution to the coronal heating and the ac- 
celeration i _oJ_th^_soJ^x_wii^ has been worked out by several 
autho r s (IDe Moortel et all 120001: iMcEwan fc De Moortell 
120061 : IDeforest fc Gurmanl Il99a 
Propagating slow waves can also be 
coronal seismology. Some authors 



O'Shea et al.l 



20071) . 



used 
used 
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properties to derive coronal paramet ers that are difficul t 
to obtain from direct observations ([Wang et al.l l2009bl ). 
But the recent detailed analyses using the spectral line 
profiles from these structures reveal periodically enhanced 
emission in blue wings, which is indicative of quasi-periodi c 
upflows (|De Pontieu fc Mcintosh! 120101 : iTian et al.l l201lh . 
It is also argued that a quasi-static plasma component in 
the line of sight may induc e the slow wave s to p roduce 
similar spectral signatures ([Verwichte et al.l l2010f ). This 
makes the g eneral nature of these dis turbances somewhat 
ambiguous. Nis hizuka fc Haral (|201ll ) reported a case of 
slow-mode waves propagating on top of a continuous 
outflow. They found that the wave signatures become 
dominant at larger heights. It is important to establish the 
nature of these disturbances before proceeding to evaluate 
their importance in heating or to other applications. To 
understand the nature of these disturbances, we study 
the distribution of power at different frequency bins for 
different open coronal structures. 
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2. Observations 

A three-hour-long imaging sequen ce taken by the 
Atmospheric Imaging Assembly (AIA; iLemen et al.ll2Q12l) 
on-board SDO, on June 5, 2010 starting from 10:00:11 UT, 
was used in this analysis. We chose two different coronal 
channels centred at 171 A, and 193 A, with their peak tem- 
perature responses at 0.8 MK, and 1.25 MK, respectively. 
Images at level 1.0 were already subjected to basic correc- 
tions that involve removal of dark current, de-spiking, flat- 
fielding, and bad-pixel removal. These were then processed 
to co-align all images in different channels to a common 
centre, to correct the roll angles, and to rescale the images 
to a common plate scale using the aia_prep.pro routine 
(version 4.10), available in SolarSoftWare (SSW). The final 
pixel scale is w 0.6" and the cadence of the data is 12 s. 

Three different open-loop structures, one isolated on- 
disk plume-like structure, a set of fan loops extending off- 
limb, and plume/interplume structures above the north 
pole of the sun were chosen in this study. Square field sub 
regions of dimensions 350, 600, and 700 pixel units, respec- 
tively, were considered to cover these three target regions. 
In both channels, images corresponding to the on-disk re- 
gion were tracked for solar rotation and then co-aligned us- 
ing the intensity cross-correlation, taking the first image as 
reference. Images in the off-limb locations were co-aligned 
directly because it was not possible to correct for solar ro- 
tation. 

3. Analysis and results 

We performed a wavelet analysis for the time series at 
each pixel to determine the distribution of power at dif- 
ferent frequencies. Keeping the total length of the time 
series (3 hrs) in mind, we decided to eliminate periods 
longer than 30 min. This was accomplished by defend- 
ing the time series with a 150 point (« 30 min) smoothed 
original light curve before performing the wavelet analy- 
sis. Frequencies with power above the 99% confidence level 
were chosen to be signif icant for a white noise process 
([Torrence fc Compolll998f ). We then considered three dif- 
ferent periodicity ranges, which we refer to as short (2 min 
- 5 min), intermediate (5 min - 12 min), and long (12 min - 
25 min), and produced power maps corresponding to these 
by adding the total significant power within the periodicity 
limits at each pixel location. Figure [T] displays the original 
intensity image (top left panels) along with powermaps for 
these periodicity ranges in the 171 A (top) and 193 A (bot- 
tom) channels of AIA for the off-limb loop structure. The 
powermaps clearly show that the power in longer periods 
is significant up to larger distances along the loop. In other 
words, this could mean that the longer periods travel far- 
ther compared to the shorter ones, which become damped 
faster. This trend is followed in both channels. This is true 
even for the other two structures, the on-disk plume-like 
structure and the north polar plume/interplume structures. 
The powermaps for these structures are shown in Fig. [5J 
Since this behaviour is the same, it is possible that the na- 
ture/source of these disturbances is the same. 

We now focus on the on-disk plume structure for a more 
detailed analysis. We traced a portion of the loop, con- 
structed space-time maps and processed them to enhance 
the visibility of the alternating ridges. Fig. [3] shows the lo- 
cation of the loop on a snapshot displaying the plume like 
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Fig. 1. Powermaps in different periodicity ranges for the 
off-limb loop structure in the 171 A (top) and 193 A (bot- 
tom) channels of AIA. The intensity image is shown in the 
top left panels. All powermaps show the total power above 
the 99% confidence level in the periodicity limits. 



structure along with the processed space-time maps in the 
two AIA channels. The processing was performed by sub- 
tracting a background trend constructed from a 55 point (^ 
11 min) smooth average of the original and then dividing 
by it. The trend subtraction removes or suppresses periods 
longer than 11 min, and division by the trend compensates 
for the intensity fall along the loop. We measured the prop- 
agation spe ed at three clean ridges follo wing the method 
described in lKrishna Prasad et al.l (|2012f ). This method in- 
volves identifying the peak positions in time from the space- 
time map at each spatial location along the ridge and then 
fitting them linearly. The peak positions were identified 
from the locations of maxima within a strip wide enough 
to cover the ridge alone. The apparent propagation speed 
and the error in its estimation could then be derived from 
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Fig. 2. Top row: Powermaps constructed in different periodicity ranges for the on-disk plume structure in 171 A (left) 
and 193 A (right) channels of AIA. The intensity image is shown in the top left panels. All powermaps show the total 
significant power above the 99% confidence level in the periodicity limits. Bottom row: Same, but for the north polar 
region 



Table 1. Periodicity and propagation speeds in different 
AIA channels. Second-strongest periods are also listed in 
brackets. 



Channel 


Periodicity 


Apparent speed (km s i ) 


(A) 


(min) 


Ridgel Ridge2 Ridge3 


171 


9.4 (15.7) 


88.7±9.7 80.5±5.7 70.6±5.8 


193 


9.4 (15.7) 


100.8±6.4 84.0±8.6 85.0±13.3 



the slope of the fit and the goodness of the fit. The bound- 
aries of the chosen strips and the final fitted peak positions 
are shown as dotted and solid lines along the ridge, respec- 
tively, in the space-time maps shown in Fig. [3) Unlike the 
case of choosing two points manually on each ridge, this 
method takes several points along the ridge and uses the 



peak positions in the spee d estimation and he nce improves 
the error in it. Recently, iKiddie et al.l ([20121 ) used a sim- 
ilar method to determine the speeds of such propagating 
disturbances and discussed their temperature dependence. 
lYuan fc Nakariakovl (|2Q12f ) also developed three different 
methods for measuring the apparent speed of propagating 
waves with considerably high accuracy. Apparent propaga- 
tion speeds obtained from the two AIA channels are listed 
in Table [U The same ridges are chosen in both channels 
for comparison. The observed apparent propagation speed 
is consistently higher in all three ridges in the hotter 193 A 
channel. To estimate periodicity, we performed a wavelet 
analysis at an arbitrarily chosen row marked by a horizon- 
tal dashed line in the space- time maps. The wavelet plots at 
this location in the 171 A and 193 A channels are shown in 
Fig. |U The first two significant periods obtained from this 
analysis are labelled at the middle right portion of these 
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Fig. 3. Left: Snapshot displaying the on-disk plume-like structure. The white lines overplotted bound the chosen loop 
location. Right: Processed space-time maps constructed from this loop portion in the 171 A (top) and 193 A (bottom) 
channels. The overplotted slanted solid lines connect the final peak positions derived from the linear fit. The slopes of 
these lines give the propagation speeds. The horizontal dashed lines in each panel mark the location of the row used in 
the wavelet analysis for periodicity estimation. 
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Fig. 4. Wavelet plots for the time series from the row marked by the dashed line in Fig [3] in 171 A (left) and 193 A 
channels. In each of these figures top panels show the original light curve in solid line and the subtracted trend in dashed 
line. The middle panels display the trend subtracted and normalized light curves. The main wavelet plot is shown in the 
bottom left and the global wavelet plot in the bottom right panels. Overplotted dotted line in the global wavelet plot is 
the 99% significance curve. Top two strongest significant periods are listed as PI and P2 in the middle right portion. 



plots. These periodicity values are also listed in Table [TJ 
As can be seen from these values the periodicity of these 
disturbances is 9.4 min. The values listed in brackets in the 
table are the second-strongest peaks. Note the presence of 
the longer periods despite the 11 min filtering. They are 
only suppressed but not completely eliminated. 

We also investigated the spatial damping along the loop. 
Fig. [5] shows the intensity values along the loop (solid grey 
line) relative to the background, averaged over three suc- 
cessive time frames. Overplotted plus symbols in blue are 
the 10 point (« 4.3 Mm) smooth average values. A damping 



sine wave function given by 

1(1) = A sm(ujl + 0)exp(-Z/L d ) + B + B ± l (1) 

was fitted to these values, which is also overplotted as a 
solid line in red. Here Aq is the amplitude of the sine 
wave, I is the length along the loop, lj and (j) are frequency 
and phase of the wave, and L^ represents the damping 
length. The last two terms were used to fit the background 
trend, if present. The first-order fit is found to be more 
than sufficient for this purpose. The best fits were obtained 
using a Levenberg-Marquardt least-squares minimisation 
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Fig. 5. Top: Solid line in grey represents the intensity val- 
ues along the loop, relative to the background, averaged 
over three successive time frames in 171 A channel. Plus 
symbols in blue are the 10 point (~ 4.3 Mm) smooth aver- 
ages of these values. The damping sine curve fitted to these 
smoothed values is overplotted as a red solid line. Obtained 
fit parameters are labelled. Bottom: Similar curves and val- 
ues generated from the 193 A channel. 



(|Bevingtonlll969f ) . Top and bottom panels of Fig. [5] display 
the plots from the 171 A and 193 A channels, respectively. 
Obtained fit parameters are labelled in the respective pan- 
els. Clearly, the relative amplitude of the oscillation and 
the damping length (L^ value in equation [T]) are lower and 
the wavelength is slightly longer in the hotter 193 A chan- 
nel. The damping lengths are 38.9 Mm and 24.5 Mm in 
the 171 A and 193 A channels respectively. Lower damping 
lengths indicate faster damping in the 193 A channel. 



4. Slow-wave model 

The farther propagation of the lower frequencies and the 
stronger dissipation in hotter channels may be understood 
as an effect of the damping of slow magneto-sonic waves 
by thermal conduction. We consider a sim ple ID model 
(|De Moortel fc Ho od 2003: iOwen et al.ll2009l ). where radial 
structuring has been n eglected. M o re co mplicated models 
have been studied by ISoler et all ([20081 ). The dispersion 
relation of slow waves in such a simple configuration is given 
by 

uo 3 - ijdc 2 k 2 uo 2 - ujk 2 c 2 + idk 4 cj = 0, (2) 



where uj is the angular wave frequency, k is the wave num- 
ber, 7 is the ratio of specific heats and c s is the sound speed. 
The parameter d = ^ 7 1 ) m: ii To expresses the importance of 
the thermal conduction. It has the dimension of time. The 
quantities po and To are the background pressure and tem- 
perature, respectively, ku = tzoT^ 2 is the thermal conduc- 
tion coefficient. Low d means a weak thermal conduction, 
high d corresponds to strong thermal conduction. Since the 
observed waves are presumably driven from near the foot 
points of the loop, we considered the frequency uj to be 
fixed by the driver, and solved the dispersion equation for 
the wave number k. For the vanishing thermal conduction 
(duj = 0), we find k 2 = uj : 



sound wave solution. In the limit of weak thermal conduc- 
tion (duo <C 1), the solution may be found as 

1 



uj f duj . \ uj 



^cond 



(3) 



where we defined the thermal conduction damping length 
^cond- This formula is only marginally applicable in the so- 
lar corona, si nce the thermal conduction in coronal loops is 
quite strong dVan Doorsselaere et al.ll201ll). For instance, 



! /Cg, which corresponds to the 



the results from IVan Doorsselaere et al.l ([201 if ) indicate 
duo w 0.7 which is less than 1 but not very low. However, 
a lot of the physics relevant for these observations may be 
understood from this asymptotic formula. 

We first consider the effect of temperature. When the 
temperature increases (i.e. when we observe the structure 
in a higher temperature passband), the thermal conduc- 
tion coefficient k\\ increases as well. This results in a higher 
value for and thus a lower value for the damping length 
L con d in the limiting case of weak thermal conduction 
(duo <C 1) that we currently consider. As found in the 
observations, the slow waves are more heavily damped in 
higher temperature passbands. Now, we consider the fre- 
quency dependence of the damping. In Eq. [3j we have de- 
fined L con d = 2c s /duo 2 ( / y — 1), where it is now obvious 
that the damping length is inversely proportional to the 
frequency squared, or directly proportional to the period 
squared. Indeed, in the observations we find that the longer 
periods propagate much higher (longer damping lengths). 
This simple model explains the qualitative features of the 
observations quite well. This model could be tested by mea- 
suring the damping lengths for different frequency bands, 
and showing that a power law with slope 2 is found. In 
principle, it also possible to use this simple model to mea- 
sure the thermal conduction. By measuring the damping 
length in a loop for several frequency bands, the thermal 
conduction may be found as the intercept of the power law 
between the frequency and damping lengths. 

5. Conclusions 

Open-loop or plume-like structures both off-limb and on- 
disk show that propagating disturbances with longer peri- 
ods travel farther distances in the powermaps constructed 
in different periodicity ranges. The apparent speeds ob- 
tained for the propagating disturbance in the on-disk plume 
structure are higher in the hotter (193 A) channel. The am- 
plitude of the disturbance is relatively low in this channel. 
The disturbance is also observed to be damped along the 
loop. The damping lengths were found to be 38.9 Mm and 
24.5 Mm in 171 A and 193 A channels, respectively. The 
faster damping of the disturbance in hotter channels might 
indicate possible damping due to thermal conduction and 
matches the earlier theories on slow-mode wave damping. 
All these observed properties can be explained using a sim- 
ple slow-wave model. We conclude that these disturbances 
are more likely to be signatures of slow mode waves than 
quasi-periodic high-speed upflows. Furthermore, these in- 
tensity waves are omnipresent at longer periods, higher up 
even in active regions. 
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